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qq Abstract. We calculate the th order in opacity number distribution of massive 

I photons (gluons) for heavy quark production radiation including interference 

from the away-side jet. While consistent with the soft photon (gluon) approxi- 
pH mation, we find that taking 1 — x ~ 1, as done in previous calculations, strongly 

affects the magnitude of energy loss. Restoring gauge invariance by including 
the radiation associated with the away-side jet fills in the "dead cone," but is a 
relatively small effect. The Ter-Mikayelian reduction from vacuum energy loss 
is 10-40% for 5-25 GeV charm and bottom quarks. 
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1. Introduction 
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The confidence in the application of perturbative QCD methods to jet energy loss 
in heavy ion collisions gained from the early quantitative understanding of 7r and rj 
suppression with null direct-7 control [EH2] has recently come into serious doubt [[3]. 
Evidence from measurements such as high-p^ correlations [ HI [5] and nonphotonic 



X 

electrons [ El [7] demonstrates clear disagreement with perturbative models [ 
O EU HU [T5]. Several papers postulate alternative nonperturbative energy loss 
mechanisms [ 1131 1141 1151 116] . and a new measurement, the double ratio of charm 
to bottom nuclear modification factors, has been suggested as a robust observable 
for testing some of these novel ideas [ I17j. Electromagnetic radiation from quark 
jets, as it is transparent to the partonic medium, holds enormous promise as a new 
tool for investigating jet energy loss mechanisms; naively one expects the spectrum 
of photons emitted from a jet that underwent a smooth, exponential slowdown will 
differ greatly from one that suffered the emission of a few hard gluons will differ from 
one that experienced a large number of soft scatterings. Nascent data of 7-hadron 
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correlations from p + p collisions additionally motivates the theoretical exploration 
of photon bremsstrahlung in heavy ion collisions [ 18 . 

In this paper we calculate as a warmup problem, and ultimately as an interesting 
problem in its own right, the th order in opacity energy loss of a heavy quark jet, 
the radiation associated with the production of a hard parton. We will generalize 
the problem of a zero mass quark emitting a zero mass photon to a massive quark 
emitting a massive photon. The lack of theoretical consistency in the understanding 
of light and heavy flavor jet suppression makes massive quark calculations of especial 
interest [|5J 1101 111] ; moreover heavy quark predictions from pQCD will be necessary 
for comparison to AdS/CFT heavy quark drag results [HZ]. For the case of th order 
emission QCD and QED are identical but for the replacement of olem with a s and a 
color Casimir. Using a massive photon will allow a comparison to already published 
results on the QCD Ter-Mikayelian effect [HH], whose main results were: (1) the 
Ter-Mikayelian effect leads to a large reduction in th order energy loss (~ 30% for 
charm quarks) , (2) the full 1-loop HTL gluon propagator can be well approximated 
by using a fixed gluon mass m g — — /x/v2, and (3) the small- ir, soft gluon, 
number distribution for th order in opacity is 

d_N%D = Q 2 a k 2 

<Pk tt 2 w [k 2 + m 2 + x 2 M 2 } 2 ' U 

where, as usual, a bold variable represents a transverse two-vector. 



2. Calculation 



As a first step to compare to previously published results [ I19j we found the number 
distribution of emitted photons when simply plugging in massive 4-vectors into a 
standard classical E&M calculation: 



E 



EM 



(2tt) 3 2 



E 

A=l,2 



ex{k) 



P 



2p-p' 



P 



k ■ p' k ■ p 



M 2 



M 2 Q 



dN {0) 

d 3 k 



(2tt) 3 2 Uk -p'){k -p) (k-p') 2 (k-p) 2 l 
rv „ k 2 + (1 — x) 2 m 2 



[k 2 + (l-a;) 2 TO 2 +a; 2 M 2 l 2 



(2) 
(3) 
(4) 



where to get from the first to the second line we used the completeness relation 
Sa=i 2 e \ e \* ~^ ~9^v, and where we took 
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k = [xE + , 



(7) 



with brackets indicating lightcone coordinates. 

There are two elements seen in Eq. Q and not in Eq. ([IJ: several factors of 
(1 — x) 2 , and an m 2 in the numerator. The first makes no difference in the limit of 
small x; however energy loss calculations integrate over all x, and it turns out that 
neglecting these factors is a large effect. The second simply cannot be reconciled 
with Eq. (|T|) . Interestingly this extra mass term in the numerator fills in the "dead 
cone," the region of small angles with respect to the jet axis for which dN g /dx — > 
in Eq. ([I]) as k — > when Mq ^ 0; this motivates additional study because naively 
the dead cone leads to a reduction in heavy quark energy loss, inconsistent with the 
observation of similar suppression patterns for pions, decay fragments from gluons 
and light quarks, and nonphotonic electrons, decay fragments from heavy charm 
and bottom quarks. 

One may rightly object that the results of Eq. Q were derived using the usual 
massless photon E&M formulae. Surprisingly the only modification of Eq. Q when 
using the Proca Lagrangian is to change the polarization sum to include the longi- 
tudinal mode. It turns out that the extra terms generated by the application of the 



identity £a=i,2,3 e A e A* = ~9\u> + k^K/r 
for massive photon radiation. 
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exactly cancel, and Eq. (W| is also valid 
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Fig. 1. The two diagrams contributing to the th order in opacity photon/gluon radiation 
spectrum. Note the inclusion of the radiation from the away side jet, which is usually 
ignored in pQCD calculations. 



In order to understand the discrepancy from the field theory perspective, con- 
sider the diagrams contributing to the th order shown in Fig. [l] Evaluation of 
these leads to 



iM = Qeu(p) 
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v(p'), 



(8) 



where we have taken A4o(p + k,p') w Aio(p,p' + k) sa Ai in the small x, soft 
radiation limit. Most pQCD calculations ignore the away-side jet; one can easily 
see that the second term in Eq. (pi), corresponding to the inclusion of the second 
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diagram in Fig.[T] is crucial for preserving the Ward identity [EQ]- Simultaneously 
dropping the ft in the numerator and the m 7 in the denominator (consistent with 
the soft photon limit) exactly reproduces the classical Proca result. We note that 
assuming Mq commutes with £*ft and retaining m 7 =/= in the denominator of 
Eq. (JsJ results in a dN p Qcv/d 3 k with leading order identical to Eq. Q but with 
(1 — x) 2 — > (1 — x/2) 2 as the prefactor of the m 7 in the numerator and (1 — x) 2 — > 
(1 — n) as the prefactor of the m 7 in the denominator. 

3. Size of Effects 

We wish to investigate quantitatively the effect of these extra terms on the th order 
energy loss. To do so we enforce physicality by restricting the x and k integration 
limits so that the emitted photon has E 1 > m g and leaves the jet with Ej et > Mq. 
For ease of comparison with [fTJ] we set /i = .5 GeV and a = .5 fixed. One can see 
from Fig. [2] the large (50-150%) effect on AE/E of including the overall prefactor 
of (1 — x). Filling in the "dead cone" makes only a small difference to the energy 
lost (5-20%); this is a surprise as the "dead cone" is the usual naive justification for 
heavy quarks having smaller radiative energy loss than light quarks. 




(a) (b) 

Fig. 2. (Color online) th order radiative energy loss for (a) charm and (b) bottom 
quarks. All results are to leader order (LO) in 1/E + . One sees that the largest effect 
(50-150%) comes from including the (1 — x) 2 prefactor and that filling in the "dead cone" 
with the massive photon is a rather small one (5-20%). Comparison with m 7 = yields 
the magnitude of the LO Ter-Mikayelian effect (10-40%). 

Fig. [3] demonstrates the effect of including all the terms generated by Eq. Q, 
not just the LO in 1/E + . Of course at higher E and pt the additional terms make 
little difference, but they regulate the otherwise divergent results in Apx/pr as 
Pt — > 0. The Ter-Mikayelian effect, given by the difference between the m 7 7^ 
and m 7 = plots in Fig. [3l varies from 10-40% for charm and bottom energy loss. 
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Fig. 3. The effect of including all terms from Eq. Q instead of just the leading order (LO) 
terms in 1/E + for (a) AE/E and (b) Apr/pr (the legend in (a) applies to both plots). For 
AE/E, the size of the relative difference in magnitude-the Ter-Mikayelian effect-is changed 
little while the overall normalization is significantly altered at low energies. For Apr/pr 
both the relative and overall normalizations change quite a bit, with the inclusion of all 
terms regulating the pr — * divergences in the vacuum production radiation spectrum. 

4. Conclusions 

Unfortunately after many years of effort there is still no single satisfactory energy 
loss model for heavy ion collisions at RHIC. This leads to the need for basic experi- 
mental tests of the gross features of the underlying energy loss mechanism, whether 
it be more like pQCD, AdS/CFT, or some other approximation. Medium induced 
photon bremsstrahlung has the potential to provide unprecedented insight into the 
modes of jet energy loss, and in this paper we took the intermediate step of ana- 
lyzing the th order in opacity production radiation energy loss. While enforcing 
gauge invariance by not neglecting the away side jet fills in the "dead cone," this 
ultimately has only a small effect on the radiation spectrum. On the other hand 
neglecting the overall factor of (1 — x) 2 in the emitted photon distribution makes 
a surprisingly large difference. This prefactor, also neglected in medium-induced 
gluon radiation derivations [|5TJ [55], may significantly alter Raa(j>t) calculations, 
especially at smaller momenta. 
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